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Abstract 

In this research article, we discuss the non-steady state concentration of mediator of electrochemical induction in chemical 
reactions. This paper presents an approximate analytical method (Homotopy and new approach to the Homotopy 
perturbation method) to solve the nonlinear differential equations. A simple and closed form of an analytical expression of non- 
steady state concentrations and effectiveness factors of current is pertaining for all values of the dimensionless parameters. The 
numerical simulation of the problem is also reported using the Matlab/Scilab software program. The analytical results are 
compared with numerical results and are found to be good in the agreement. The obtained results are valid for the whole 
solution domain. 

Keywords 

Mediated Concentration; Current; Mathematical Modeling; Initial-boundary Value Problem; Numerical Simulation; Laplace-homotopy 
Perturbation Method 

Introduction 

There are many methods for analysis electro-chemical reactions. One of the methods is a mediated electro analysis 
by use of electrochemical induction, that is used in this paper. The redox compound mediated and electro 
catalyzed methods process is called mediated electro analysis [Medebille et.al (1991)]. Mediated electrochemical 
induction is indirect and also more systematically exploited [Anderiex et. A1 (1980)]. The main idea of this problem 
is electron transfer in reactant to product with mediators [Medebille et.al (1997)]. In addition redox titration 
between the substrate and the result of the experiment. The electrochemical potential of redox active compounds 
can be determined in solution by applying a redox titration.The absorption spectrum of the investigated compound 
is observed [Medebille et.al (1991)]. The potential of the solution is changed by adding a reducing or oxidizing 
agent. If the potential in the solution is negative (positive) enough to reduce (oxidization) the investigated 
compound, the optical spectrum will change until the reduction (oxidation) is complete [Andrieux et.al (1991)]. In 
the mediated electrochemical process is faster than in, direct electrochemical induction [2], In this paper the author 
analysied electrochemical induction of the nucleophillic substitution of perfluorohexyl iodide by 2-methyl-5- 
nitroimidazolate ions (Nu ) in acetonitrile mediated by the 4-nitropyridine-N-oxide anion radical, and give an 
experimental by use of above described method [Medebille et.al (1991)]. The indirect (mediated) electrochemistry 
should be more systematically exploited. Areactant is generated electrochemically and its homogeneous reaction 
with the substrate can be monitored by the ensuing variation of the current . The reactions of single electrons as 
well as of electron pair donors or acceptors can be followed by this method [Bard et.al (2001)]. This in 
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situgeneration allows the investigation of very reactive species, the only limit being that they should react with the 
substrate faster than they react with any other molecule present in the solution [Savant et.al (1997)]. Then 
mathematical modeling eqns (1-5), are converted in theose eqns. In the mediated electrochemical process these 
reactions are faster than electron transfer to (or) from B , the overall kinetics will depend on a competition between 
reaction and first order deactivation reaction [Medebille et.al (1991)]. The purpose of this paper is to derive the 
analytical expressions for mediated concentrations for substrate and the dimensionless current using the new 
approach to Homotopy perturbation method. 

Mathematical Formulations of the Problem 

The mathematical model for estimating the electrochemical reactions was set forth by Saveants (1997) and further 
developed by Medebille et. al (1991). This problem based on electrochemical reaction with the use of redox titration. 
The substrate reactant converted to products using redox couple P/Q [Medebille et.al (1991)]. 
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(3) 
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The above kinetics reaction eqns.(l)-(5) are formed by the following set of nonlinear partial differential 
equations, with initial and boundary conditions: 


dp dp 

— = — +\aq-X i cp 
or oy~ 

(6) 

dq d 2 q 

- = — -\aq + l 3 cp 
ox dy~ 

(7) 

da d 2 a 

(8) 


where 


X 2 b = X 3 cp, \aq = (1, + X 4 )b 


The initial and boundary conditions are, 

r = 0, y > 0 and y = <x> z>0:p = l,q = 0,a = y 


(9) 

( 10 ) 


„ ^ „ op Gq n , da _ 

y = 0,x>0:— + — = 0,p = q*e , — = 0 
dydy dy 


( 11 ) 


where 


t = (FV/RT)t, £ = -(F/RT)(E-E° pq ), 

i 

y = x(FV / RT) 2 


( 12 ) 


8 



Advances in Chemical Science (SEPACS) Volume 4, 2015 


www.seipub.org/sepacs 


4 = 
4 = 


RT k/ p 
F v 
RT k/‘ 
F v 


RT Kc° 

i — 2 P 2 — 

, /l 2 ,/l 3 

F v 



4 

4+4 


RT k,c° p 
F v 


k 2 + k 4 


(13) 


where a, p, q are the concentrations of the corresponding substrates, that's been normalized, with respect to the 
concentration of substrate of the bulk solution c ° p and k 1 , k 2 4, , k 4 are rate constants of the corresponding 
reactions. 

We introduce the a* = cr a into the precedingeqns, we obtain with the initial and boundary conditions. 


8q 

dx 



(14) 


da 

8x 


8- a , , 

— ; V 4 

d-y 


(15) 


r = 0, y > 0, y = oo, z > 0 : q = 0, a =y =ay 


(16) 


y = 0, r > 0 : q 


-L-, ^ = o 

l + e' dy 


(17) 


The cyclic voltametric wave is given by. 


= (18) 
I\ 0.446 [dy)^ 

Analytical Solutions of the Concentrations and Current Using the Homotopy and New 
Approach to Homotopy Perturbation Method 

Recently, many authors have applied the Homotopy perturbation method (HPM) to solve the non-linear problem 
in physics and engineering sciences [He (1999 & 2003), Ariel (2010), Ananthaswamy et.al (2013) and ]. This method 
is also used to solve some of the non-linear problem in physical sciences [Shanthi et.al (2013 & 2014) ]. This method 
is a combination of Homotopy in topology and classic perturbation techniques. Ji-Huan He used to solve the 
lighthill equation, the duffing equation and the blasius equation (1999 and 2003)]. The HPM is unique in its 
applicability, accuracy and efficiency. The HPM uses the imbedding parameter p as a small parameter, and only a 
few iterations are needed to search for an asymptotic solution. By solving the eqns. (14) and (15) the mediated 
concentrations q(y,t) and a(y,t) can be obtained by the following eqns. 


q(y,t) 


erf 


-tslb 


2 tyfb _y 

2 4t 
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+ e 2y ' rb erfc 


2 t\fb + y 

2 V7 


2(e~ f + 1) 


(19) 
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Vn 

We use the value of mediated concentration <7 in the eqn.(20), then we get the normalized current inthe eqn.(21). 

vA 


r dq N 


<9y 


V'UVjM, 
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( 21 ) 


- Vd etf (l4bt) + - l^b erfc(2.4b + 1) 

4t(e~ 4 +1) 


Numerical Solution 

The nonlinear differential eqns. (14) - (17) are also solved using numerical methods. The function pdex4 in 
Matlab/Scilab software function which is solving the partial differential equation. Its numerical solution is 
compared with Homotopy perturbation method and it gives satisfactory results. 

Results and Discussion 

Figure 1 and 2 show the mediator concentration q( y,t ) versus dimensionless distance y. From Fig. 1 (a), it is clear 
that when the the excess parameter y increases, the corresponding mediator concentration q(y,t) decreases in 
some fixed values of the other dimensionless parameters cr, /, , , c and t. From Fig. 1(b), we notice that when the 
ratio of rate constant <3 increases, the corresponding mediator concentration q( y,t ) also increases in some fixed 
values of the other dimensionless parameters. From Fig. 2(a), we infer that when the spherical diffusion parameter 
E increases, there is no change in the mediator concentration q( y,t) in some fixed values of the other dimensionless 
parameters. From Fig 2(b) it is clear that when dimensionless time t increases, the corresponding mediator 
concentration q( y,t ) also increases in some fixed values of the other dimensionless parameters. From Fig 2(c) is 
observe that when the rate constant decreases, the mediated concentration q( y,t ) decreases in some fixed 
values of the other dimensionless parameters. 

Figure 3 and 4 show the mediator concentration a( y,t ) versus dimensionless distance y, for various dimensionless 
factors o, y. A, , c, t. From Fig. 3 (a), it is clear that when the the dimensionless time f decreases the corresponding 
mediator concentration a( y,t ) increases in some fixed values of the other dimensionless parameters. From Fig. 1(b), 
we notice that when the excess parameter y increases, the corresponding mediator concentration a( y,l ) increases 
in some fixed values of the other dimensionless parameters. From Fig. 4(a), its clear that when the rate 
parameter Xj decreases, the corresponding mediator concentration a( y, t ) increases in some fixed values of the 
other dimensionless parameters. From Fig 4(b) it is notice that when the excess parameter y increases, the 
corresponding mediator concentration a( y,l ) increases in some fixed values of the other dimensionless parameters. 
From Fig 4(c) when the spherical diffusion parameter E increases, the corresponding mediator concentration 
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a( y,t /increases in some fixed values of the other dimensionless parameters. Figure 5 and 6 indicate the three 
dimensional diagram for q( y,t ) and a( y,t ) with respect to Xj and a and for various values of the other 
dimensionless parameters. The normalized current \p versus the dimensionless excess parameter y* are presented 
in Figures. 7(a)-(d) in some fixed values , = 0.45, c, = 0. 1 , a = 3 . Figure 7(a) and (c) it is clear that when ratio of rare 
constant , spherical diffusion parameter and dimensionless time,corresponding the normalized current decreases 
for various values of dimensionless parameters. From Fig 7(d) it is clear that the when rate parameter /, , increases 
the corresponding normalized current increases for different values of dimensionless parameters. Table 1 
represents the error table for the analytical results (eqns. (19) and (20)) and numerical simulation for the 
dimensionless parameters h, . 




FIGURE 1. PLOT OF THE MEDIATED CONCENTRATION q( y,t ) VERSUS DIMENSIONLESS DISTANCE Y IN SOME FIXED VALUES OF 
THE DIMENSIONLESS PARAMETERS y, and a WHEN (a) ^ = 0.45, % = 0.001, r = 1. a = 3 AND y = 5,7,9 

(b)Aj =0.45, ^ = 0.001, r = l,y = 4 AND a = 6,9,12 USING THE EQN. (19). 




Distance y 



Distance y 

FIGURE 2. PLOT OF THE MEDIATED CONCENTRATION q( y,t ) VERSUS DISTANCE Y IN SOME FIXED VALUES OF THE 

DIMENSIONLESS PARAMETER A|,^,f,y, and a, WHEN 
(a) A.! = 0.45, t = 1, y = 4, a = 3 and £, = 0.01. 0.05, 0.007 (b) £, = 0.001, f = 1, y = 4, a = 3 AND Aj = 5,3,1 (c) 

Aj =0.45, ^ = 0.001, y = 4,a = 3 and t = 0.05.0.01,0.007 USING THE EQN.(19). 
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FIGURE 3. PLOT OF THE MEDIATED CONCENTRATION a(y,t) VERSUS DIMENSIONLESS DISTANCE Y IN SOME FIXED VALUES OF 
THE DIMENSIONLESS PARAMETERS X^t, y, and a WHEN (a) ^ =0.45, 5 = 1, y = l,a = l AND t = 0.57,0.65,1 (b) =0.45, 5 = 1, f = l,y = l 

AND a = 0.5,1, 1.5 USING THE EQN. (20). 



Distance y 

FIGRE 4. PLOT OF THE MEDIATED CONCENTRATION a(y,t) VERSUS DIMENSIONLESS DISTANCE Y IN SOME FIXED VALUES OF 
THE DIMENSIONLESS PARAMETERS X h ^,t,y, and a WHEN (a) f = 1, 5 = 1, y = l,a = 1 AND = 0.01, 0.25, 0.45 

(b) Lj =0.45, 5=1, t = 1,0 = 1 AND y = l, 1.5, 2 
(c) Xi = 0.45, t = 1,CT = 1 and 5 = 0.005, 0.05, 0.5 USING THE EQN. (20). 

(a) 



1.5 

0.5 

Rate parameter > i 


Distance y 



FIGURE 5. THE MEDIATED CONCENTRATION OF q( y,t ) IS PLOTTED. THE GRAPH WAS CONSTRUCTED USING EQN. (19) FOR THE 

DIMENSIONLESS PARAMETERS ^ and a. 
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0.98 

0.96 

0.94 

0.92 


Rate parameter 



FIGURE 6. THE MEDIATED CONCENTRATION OF a( y,t ) IS PLOTTED. THE GRAPH WAS CONSTRUCTED USING EQN. (20) FOR THE 

DIMENSIONLESS PARAMETERS >4 and a. 


00 (b) 





FIGURE 7. PLOT OF THE NORMALIZED CURRENT ip VERSUS y IN SOME FIXED VALUES OF THE DIMENSIONLESS PARAMETERS 
and a WHEN (a) X, = 0.45, £ = 0.001. a = 4 and t = 2, 3,4 (b) A, = 0.45. t = 1, a = 4 and £ = 0.01, 0.03. 0.05 , 

(c) A.! = 0.45, f = l, £ = 0.001 and a = 3,5,7 (d)£ = 0.001, ( = 1, o = 4 and A, = 1,3, 5 USING THE EQN.(21). 

TABLE: 1A COMPARISION OF OUR ANALYTICAL RESULT OF MEDIATOR CONCENTRATION q( y, t ) VERASUS DIMENSIONLESS 
DISTANCE Y WITH THE NUMERICAL SIMULATION FOR VARIOUS VALUES OF A , THE USING EQN.(19) 

WHEN if = 0.001, t = 1, y = 4, and a = 3 


y 

=1 

A., =3 

A.; =5 

Analytical 

Numerical 

(%) of 
deviation 

Analytical 

Numericl 

(%) of 
deviation 

Analytical 

Numerical 

(%) of 
deviation 

0 

0.50025 

0.5002 

.019 

0.50025 

0.5002 

.019 

0.50025 

0.5002 

.019 

00.2 

0.23 

0.23 

0 

0.151 

0.123 

0.03 

.1 

0.1 

0 

0.4 

0.101 

0.108 

2.84 

0.045 

0.00364. 

2.5 

0.02 

0.02 

0 

0.6 

0.045 

0.047 

4.25 

0.01 

0.01 

0 

0.0047 

0.0004 

1.75 

0.8 

0.02 

0.02 

0 

0.003 

0.002 

.5 

0.001 

0.00003 

2.3 

1 

0.02 

0.02 

0 

0.002 

0.002 

0 

0.0002 

0.000007 

3.52 


Average % of deviation 

1.15 

Average % of deviation 

0.51 

Average % of deviation 

1.26 
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TABLE: IB COMPARISION OF OUR ANALYTICAL RESULT OF MEDIATOR CONCENTRATION a(v,t) 
VERASUS DIMENSIONLESS DISTANCE Y WITH THE NUMERICAL SIMULATION FOR VARIOUS VALUES 
OF A, THE USING EQN.(20) WHEN f = 1, t = 1, y = 1, and a = 1 


y 

\ = 0.01 

\ = 0.25 

\ = 0.45 

Analytical 

Numerical 

(%) of 
deviation 

Analytical 

Numericl 

(%) of 
deviation 

Analytical 

Numerical 

(%) of 
deviation 

0 

0.9024 

0.9024 

0 

0.943 

0.943 

0 

0.988 

0.988 

0 

00.2 

0.9075 

0.9075 

0 

0.946 

0.946 

0 

0.9905 

0.9905 

0 

0.4 

0.9222 

0.9222 

0 

0.9551 

0.9551 

0 

0.9931 

0.9931 

0 

0.6 

0.9441 

0.9441 

0 

0.9441 

0.9441 

0 

0.9933 

0.9933 

0 

0.8 

0.9722 

0.9722 

0 

0.9681 

0.9681 

0 

0.9966 

0.9966 

0 

1 

1 

1 

0 

1 

1 

0 

1 

1 

0 


Average % of deviation 

0 

Average % of deviation 

0 

Average % of deviation 

0 


Conclusion 

In this paper a Homotpy and new approach to Homotopy perturbation method (HPM) is successfully applied for 
solving nonlinear partial differential equation. Approximate analytical expressions of the mediated concentrations 
of substrate in the electrochemical reaction using in induction method for all values of the dimensionless 
parameters is obtained using the HPM. The results are compared with the numerical simulation and it gives 
satisfactory agreement. This method can be easily extended to solve the other strongly nonlinear initial and 
boundary value problems in physical, chemical and biological sciences. 
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Appendix A: 

Basic Concept of the Homotopy Perturbation Method [He (1999 and 2003)] 

To explain this method, let us consider the following function: 

D 0 (u)-f(r) = 0, r e Q (A.l) 

with the boundary conditions of 

du 

BJu , — ) = 0, r g T (A.2) 

dn 

where D 0 is a general differential operator, B n is a boundary operator, fir) is a known analytical function and T is the 
boundary of the domain Q . In general, the operator D 0 can be divided into a linear part L and a non-linear part N . The 
eqn.(A.l) can therefore be written as 

L(u) + N (u) - f (r) = 0 (A. 3) 

By the Homotopy technique, we construct a Homotopy v(r, p) : Q x [0,1] — > 5? that satisfies 

H (v, p) = ( 1 - p)\L(v) - L(u 0 )] + p[D o (v) - = 0 ( A.4) 

H(v, p ) = L(v) - L(u 0 ) + pL(u 0 )+ p[N (v) - fir)] = 0 A.5) 

where p e [0, 1] is an embedding parameter, and u 0 is an initial approximation of the eqn.(A.l) that satisfies the boundary 
conditions. From the eqns. (A.4) and (A.5), we have 

Hiv,0) = Liv)-Liu Q ) = 0 (A.6) 

H (v, 1) = D o (v) - fir) = 0 (A .7) 

When p= 0, the eqns. (A.4) and (A.5) become linear equations. When p =1, they become non-linear equations. The process of 
changing p from zero to unity is that of L(v) - L(t/ 0 ) = 0 to D o iv) -fir) = 0 . We first use the embedding parameter p as a small 
parameter and assume that the solutions of the eqns. (A.4) and (A.5) can be written as a power series in p : 

v = v 0 + pv , + p 2 v 2 + ... (A.8) 

Setting p = 1 results in the approximate solution of the eqn. (A.l): 

u = lim v = i' 0 + Vj + Vt + ... (A.9) 

P->1 

This is the basic idea of the HPM. 
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Appendix B: 

Solution of the Non-Linear Eqns. (14)-(17) by Using the Homotopy Perturbation Method 


Using the Homotopy and new approach to the Homotopy perturbation method to solve the eqn (14) and (15) . In this appendix 
B, we have indicated how to determine the solution of eqn. (14-15) using the boundary conditions eqs. (16) and (17). In order to 
solve eqn. (14-15), we construct the Homotopy for the equation as follows: 


o-p) 


d q J * , , dq 

77 -V 9(y=°°)-— 

oy ot 


dq . . dq 

— - -am q 

dy 2 dt 


= 0 


(B.l) 


0 ~P) 


dji 

dy 2 


dq 

8t 


dq . . dq 

— j-\a q 

dy 2 dt 


= 0 


(B.2) 


The approximate solution of eqn. (B.2) is 
Subtracting eqn (B.3) into eqn (B.2) and arranging the coefficients of powers of p we get 


q = p° % + p l q, + p 2 q 2 + • 


0 8-q 0 8q 0 

P -—r— Vt?o=° 

dy ot 

i d 2 q i dq, , 

p ■— -j — ; — h°yq\ V % = ° 

By dt 

The initial and boundary conditions in eqs.(16)-(17) becomes, at 

T = 0, </o = 0 anc l ?,• = 0 Vi = 1,2,3... 

x = 0 a 0 = Gy and a, =0 Vi = 1,2,3., 


y = 0,,q n = q. = 0 Vi = 1.2.3.. 

l + exp(-a 


y = o. 

— 5- = 0 and 

^ = 0 

Vi = 1,2, 3 

dy 

dy 

y = co 

O 

II 

O 

q, =o 

Vi = 1,2,3.... 

y = °o 

a o = oy. 

a, = 0 

Vi = 1,2,3.. 


and 


(B.3) 

(B.4) 

(B.5) 

(B.6) 

(B.7) 

(B.8) 


In Laplace plane, eqs (B.7-B.8), become, 


y — o°, 

< 7 ,, = — and 

s(l + <?"") 

(B.9) 


— oy 


y — oo 

a o = — 



s 


y = o. 

da n 

q n = 0 and = 0 

u n,. 

(B.10) 


Solving the eqns (B.7) and (B.8) and using the boundary conditions eqns.(B.9) and (B.10),we obtain the following results: 

(~yfs+\ <ry) 

do = 


s( 1 + e'"') 
- cry 


(B.ll) 

(B.12) 
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— -\oy 
a, = 


1 2(l + e - f ) 




sYs 


By taking the Laplace transform for the above eqns. we get the solutions in the text eqns. (18) and (19). 


Appendix C: Nomenclature 


Symbols 

Meaning 

T 

Temperature 

h 

Height of the catalytic peak 

c 

Height of the cathodic peak in the absence of substrate 

D 

Diffusion coefficient 

C l 

Concentration of oxidized mediator 

c° A 

Concentration of substrate 

F 

Faraday constant 

R 

Resistance 

E 

Electrode potential 

E° 
c PQ 

Standard potential of the mediator couple 

a 

Mediated concentration of substrate 


Mediated concentration of substrate 

t 

Time 

X 

Distance from the electrode 

7 

Excess parameter 

(7 

Ratio of rate constants 


rate parameter 

Z 

Dimensionless time 

K, 

rate constants of the homogenous chemical reaction 

V 

Dimensionless current 

V 

scan rate of solvent 

Y 

Dimensionless Parameter 

f 

Dimensionless parameter of spherical diffusion 

Y 

Dimensionless distance 

PQ 

Redox couple 

P 

Mediator 

A 

Substrate 

RfX 

Substitution of the halogen of perfuroalkyl 

X 

Halogen 

Rf 

Perfuroalkyl halides 


(B.13) 
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